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We identify a new mechanism for cooperative emission of light by an ensemble of dipoles near 
a metal nanostructure supporting a surface plasmon. The cross-talk between emitters due to virtual 
plasmon exchange leads to a formation of three plasmonic super-radiant modes whose radiative 
decay rates scales with A'^, while the total radiated energy is thrice that of a single emitter. Our 
numerical simulations indicate that the plasmonic Dicke effect survives non-radiative losses in the 
metal. 



Radiation of a dipole near a metal nanostructure sup- 
porting surface plasmon (SP) is attracting renewed in- 
terest due to possible biosensing applications 1]. While 
early studies mainly focused on fluorescence of molecules 
near rough metal films recent advances in near- field 
optics and in chemical control of molecule-nanostructure 
complexes spurred a number of experiments on sin- 
gle metal nanoparticles (NP) linked to dye molecules 
0, S, H, H S Hi or semiconductor quantum dots 
Emission of a photon by a dipole-NP complex involves 
two competing processes: enhancement due to resonance 
energy transfer (RET) from an excited dipole to a SP 
p^ . and quenching due to decay into optically-inactive 
excitations in the metal These decay channels 

are characterized by radiative, F'', and non-radiative, 
r"*", decay rates, respectively, and their balance is de- 
termined by the separation, d, of the emitter from the 
metal surface [H [H. The emission is most enhanced 
at some optimal distance, and is quenched close to the 
NP surface due to the suppression of quantum efficiency, 
Q = F^/ (F'' -I- F"''), by prevalent non-radiative pro- 
cesses. Both enhancement and quenching were widely 
observed in fluorescence experiments on Au and Ag 
nanoparticles SliSJl, 0, 0- In recent single-molecule 
measurements [y, 0, S] , the distance dependence was in 
a good agreement with single-dipole-NP models [l2l. [l^. 
prompting proposals for a NP-based nanoscopic ruler [8| . 

In this Letter, we identify a novel mechanism in the 
emission of light by an ensemble of dipoles located near a 
nanostructure supporting a localized SP. A typical setup 
would involve, e.g., dye molecules 0, 0, Q or quantum 
dots Q attached to a metal NP via DNA linkers. Namely, 
we demonstrate that RET between individual dipoles 
and SP leads to a cross-talk between the emitters. As 
a result, the emission of a photon becomes a cooperative 
process involving all dipoles in the ensemble. This plas- 
monic mechanism of cooperative emission is analogous 
to the Dicke effect for N radiating dipoles in free space, 
confined within a volume with characteristic size smaller 
than the radiation wavelength A In that case, 

the cooperative emission is due to photon exchange be- 
tween the emitters that gives rise to super-radiant (SR) 
states with total angular momentun 1 and enhanced ra- 
diative decay rate ~ ATJj, where FJj is the decay rate 




FIG. f : (Color online) Radiative coupling of emitters in free 
space (a) and plasmonic coupling of emitters near a NP (b). 



of an isolated dipole. In contrast, in plasmonic systems, 
the dominant coupling mechanism between dipoles is SP 
exchange, i.e., excitation of a virtual SP in a nanostruc- 
ture by an excited dipole followed by its absorption by 
another dipole, rather than direct radiative coupling [see 
Fig. [T] . Such a SP-induced coupling leads to the forma- 
tion oi plasmonic SR states that dominate the emission of 
a photon. Importantly, because the nanostructure acts as 
a hub that couples nearby and remote dipoles with about 
equal strengths, the SP-induced cross-talk is more uni- 
form throughout the ensemble, as compared to the radi- 
ation coupling, resulting in a more efficient hybridization 
and hence cooperative emission. 

The usual photonic Dicke effect can be suppressed by 
internal non-radiative processes in molecules or by their 
energy exchange with the environment. For example, 
when an ensemble of emitters is located nearby a metal 
nanostructure, the photon exchange processes are largely 
quenched by ohmic losses in the metal. In contrast, as we 
demonstrate in this paper, the SP exchange mechanism 
is not significantly affected by quenching up to very small 
distances, and thus provides the main channel for coop- 
erative emission in plasmonic systems. Specifically, we 
show that for an ensemble of N dipoles distributed in a 
solid angle around a metal NP, there are three plasmonic 
SR states with radiative decay rates jj^ ~ iVF'"/3. Fur- 
thermore, in a wide range of dipole-NP distances, their 
non-radiative decay rates also scale as 7^'' ~ iVF"''/3, so 
that the SR quantum efficiencies are essentially the same 
as those of individual dipoles near a NP. As a result, the 
total energy radiated by an ensemble, W, is only thrice 



2 



that radiated by a single dipole near a NP, Wq: 
W ~ 3{hkc/4)Q = 3Wo, 



(1) 



where k and c are wavevector and speed of light, and the 
remaining energy is dissipated in a NP via sub-radiant 
states. The fact that, in plasmonic systems, radiated en- 
ergy of an ensemble is nearly independent on its size could 
allow unambiguous determination of single-molecule de- 
cay rates in situations when a large but uncertain number 
of molecules participate in the emission. 

Theory — We consider a system of N emitters, e.g., 
fluorescing molecules, with dipole moments = 'ij^jj 
where dj and ej are their magnitudes and orientations, 
respectively, located at positions around a spheri- 
cal NP of radius R in a dielectric medium with its 
center at origin. We assume incoherent emission, i.e., 
molecules initially excited by a laser pulse, subsequently 
relax through internal transitions before emitting a pho- 
ton, and adopt classical model of Lorentz oscillators with 
random initial phases. The frequency-dependent electric 
field, E(r,ix'), created by all dipoles in the presence of a 
NP, satisfies Maxwell's equation 



^^^'^^^' E(r,c.)-VxVxE(r,c.) = 



■Kr,^), (2) 



where dielectric permittivity e(r, uj) is that of the metal 
inside NP, e{uj), for r < R, and that of the outside di- 
electric, eo, for r > R. Here j(r,cj) — —i e^'^*j{t)dt, 
is the Laplace transform of dipole current — 
(l'}2,jdj{t)ej5{Y — Tj), where dipole displacements are 
driven by the electric field at dipoles positions. 
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with the initial conditions dj = d^ej sin ipj, dj = 
LOodoBj cos (pj, and E = for t = (dot stands for 
time-derivative). Here ujq, q, to, and ipj are oscillators 
frequency, charge, mass, and initial phase, respectively 
{ajQ = hq^/md^). Closed equations for dj{uj) are ob- 
tained by Laplace transforming Eq. ^ with the above 
initial conditions and then eliminating E from Eqs. (|2|3p 
[l7| . The latter can be expressed via normalized displace- 
ments, Vj{u!) — dj{uj)/do — i {uJa/uj'^) cosLpj — ui^^ sinipj, 
and voj = — i ((Jq/cj^) cosipj — (^ujq/u;^ s'mipj, as 



E(r,c^) 



AirdoqajQ 
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where G(r,r',(jj) is the electric field Green diadic in the 
presence of NP. For the photon frequency close to those 
of dipoles, cj ~ we arrive at the following system. 



2^0 2 



(5) 



where the complex self-energy matrix, Sjfe, is given by 
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The system ([5]) determines eigenstates of N emitters 
coupled to each other via radiation field and electronic 
excitations in a NP. For jr^ — r^l <C A, we can use 
the near-field expansion of the Mic theory Green diadic 
G(r'j,rfe;w) [l^ for calculation of Sj^. The self-energy is 
dominated by imaginary part that contains the SP reso- 
nance. The details will be given elsewhere [l^, and the 
decay matrix, Tjk = — ImEjfc, is a sum of radiative and 
nonradiative terms, Tjk = F^^ + F"^ + SjkTQ^, where 
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and Fg"" accounts for internal molecular transitions {Sjk 
is Kronecker symbol). Here Fg = 2dgk^/3heo is dipole 
radiative decay rate (fc = ^JT^iojc is wavevector), aiiiS) = 

J^^' + i[e(c^)-eo] 



£(w)-|-(l-|-l/Oeo 
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polarizabilities, and matrices T^^ and if -j] are 
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where ■0;m(r) = V [r ' ^ !!„(?)] and X/m(r) = 
V [r'y/m(f)], Yim{Y) being spherical harmonics. Natu- 
rally, only the dipole (/ = 1) term contributes to F^^,, 
while F"^ includes all angular momenta. From diagonal 
elements, single-dipole-NP rates can be easily recovered 
for normal (s =J-) and parallel (s =|j) orientations with 
respect to the NP surface Tl ~ 



Tl\l + asai/rl\ 



and TY - (3F5/2fc3)^^6«a;Vr2 
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Radiated energy in the unit frequency interval is ob- 
tained by integrating spectral intensity over solid angle, 
dW/duj — (ceo/47r^) / |E(r,a;)|^ r'^dft, and averaging the 
result over initial phases of oscillators, ipj. Here the far- 
field E(r,w) is given by Eq. ([4]), where vj is the solution 
of Eq. ^ and G{r_^,uj) is the large r asymptotics of the 
Mie Green diadic |13| • The details will be given elsewhere 
llSl. and the final result reads 
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In the absence of dipole coupling, i.e., for purely diagonal 
Sjfe = —iSjkT and F^j, = SjkT"^, the frequency integration 
recovers radiated energy of N isolated dipoles near a NP, 
W = N{^huJo/A)Q = NWo. 

To illustrate the effect of SP coupling between emit- 
ters, first consider N dipoles randomly distributed in a 
solid angle around a NP at the same distance d > R from 
its surface, with normal or parallel orientations. At such 
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distances, the high angular momenta (/ > 1) contribu- 
tions to r^l are suppressed, and decay matrices Eq. ([7]) 
take simple form FJ^, — T^Ajk and F"^ = F"'"Ajfe, where 
Ajk = Bj -Sk is cosine matrix, and F"*" includes only I = 1 
term (s =_L, ||). We now introduce cooperative decay ma- 
trices as = (iVr-/3)5^, and 7;^: = iNrT/3)B^,, 
where B^^ = (3/-/V) Cj^ejv is 3 x 3 matrix in coor- 
dinate space with TrB = 3. Now we note that, since 
Tri" = Tr(7Vi?/3)" for any integer n, the TV x iV matri- 
ces T^j'^^ have only three non-zero eigenvalues coinciding 
with those of matrices 7^^'"''. Therefore, only these eigen- 
values contribute to the spectral function, 
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where 7^ = Xf^NT^/S and 7^ = X^N (F', + F'-) /3 
are radiative and total decay rates of plasmonic SR 
states, are their quantum efficiencies, and ^ 1 
are eigenvalues of B^^. Importantly, both radiative and 
nonradiative rates of SR states are enhanced by factors 
^ N/3 (for each degree of freedom). However, these fac- 
tors effectively cancel out in the quantum efficiencies, Q^. 
Furthermore, the contribution of Fg"" is the denominator 
of Eq. pT|) is suppressed by the factor ^ 1, i.e., 
are not sensitive to intramolecular relaxation processes. 
Not too far from a NP, when FJJ'' < (F^ + F^''), SR and 
single- molecule efficiencies essentially coincide, ~ Q. 
Integrating Eq. (jTU]) over frequency, we obtain Eq. ([1]). 

The origin of three plasmonic SR states is that even 
though dipole orientations may be uniform with respect 
to the curved metal surface (e.g., normal to it), they are 
not uniform in space, and vice versa. Note that for purely 
radiative coupling and for uniform dipole orientations, 
there is only a single SR state [3, EM EB] ■ 

Numerical simulations and discussion — Although 
Eqs. (jlOllip were derived for moderate distances, d> R, 
these results apply even close to NP surface. In Fig.[2l we 
show numerical simulations of the eigenvalue distribution 
of FJj. -I- F"^ for 30 molecules with normal dipole orienta- 
tions randomly placed in spherical angle around Au NP of 
radius i? = 16 nm (in a medium with eg = 1.77). Usually, 
stretching and folding of linker molecules causes fluctu- 
ations of dipole-to-surface distances by some amount A 
around the average value d [1, 0, so 10% spread in 
distances (A/d = 0.1) was included. Calculations were 
performed at SP energy of 2.31 eV, and NP polarizabil- 
ities ai{Lu) with up to ^ = 30 in Eq. ([7]) were included 
in F"^'. To account for quantum-size effects in small NP, 
we incorporated in ai{uj) the Landau damping of Z-pole 
plasmons, characterized by rate 7; ~ Zlvp/AR, where 
fi? is the Fermi velocity in the metal. The distribution 
of decay rates reveals the increasing role of nonradiative 
processes as the average distance to NP surface is reduced 
(see Fig. [5]). For d = 20 nm, there are only three non- 
zero eigenvalues corresponding to SR states, in agreement 
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FIG. 2: (Color online) Distribution of decay rates for 30 
dipoles around Au NP at several average (with 10% fluctua- 
tions) distances to its surface. 
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FIG. 3: (Color online) Eigenvalues of (a) radiative, F^'j., and 
(b) full, Vjk, decay matrices vs. average distance to NP sur- 
face for 30 dipoles randomly distributed around Au NP. 



with Eq. pH)) . With decreasing d, the remaining — 3 
sub-radiant states start emerging (d = 10 nm) and, at 
small distances (d = 5 nm), all system eigenstates ac- 
quire a finite decay rate. Note that a similar result holds 
for any dipole orientations; e.g., for random orientations 
in tangent plane, the three SR modes in Eq. (fTU]) are well 
separated from the rest. 

The crucial distinction between the photonic and plas- 
monic Dicke effects stems from non-radiative coupling 
between emitters in the latter. In the absence of a NP, 
the system eigenstates are eigenvectors of radiative decay 
matrix, Fg, and represent super- or sub-radiant modes 
characterized by strength of their coupling to the ra- 
diation. In the presence of a NP, optically bright and 
dark states are also defined through F'", Eq. ([7]), whose 
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FIG. 4: (Color online) Combined q uantum efficiencies for 
ensembles of 30, 60, and 100 dipoles compared to that for an 
isolated dipole near Au NP. 

eigenvectors describe three SR and — 3 sub-radiant 
modes regardless of dipole-surface separation, as shown 
in Fig. [3lja). However, the true system eigenstates, 
are described by the full decay matrix, 1"" + F'"", with 
eigenvalues Tj , whose nondiagonal elements include non- 
radiative coupling. Therefore, the true radiative decay 
rates are given by expectation values FJ = (j|F''|j), 
while the quantum efficiencies are Qj — T^/Tj. In the 
case, e.g., of normal or tangential orientations, the dipole 
terms in the non-radiative decay matrix Eq. ([7|) possess 
the same symmetry as the radiative decay matrix, so for 
distances not too close to the NP, the SR modes are still 
system eigenstates. At the same time, the high-/ terms in 
F"'" have different symmetry than F*" and, therefore, they 
cause mixing of SR and sub-radiant modes [see Fig.[3]Jb)]. 
However, except for very small distances, this mixing is 
weak so that the emission remains cooperative. In Fig.|4l 
we compare distance dependence of combined quantum 
efficiencies, Qens = J2j Qj^ ^O'' ensembles of 30, 60, and 
100 molecules to the single-molecule Q near a i? = 16 nm 



gold NP. For distances d > R/2 (8 nm), all ensemble de- 
pendences collapse into single curve with amplitude 3Q 
(we used FJf = 1.08 x 10^ s^^ for Cy5), indicating that 
the emission is dominated by SR modes. Even closer to 
NP surface, up to d « 5 nm, the emission remains co- 
operative, although deviations from 3Q behavior appear. 
For smaller d, the eigenstates are no longer SR and sub- 
radiant modes, and cooperative emission is destroyed by 
non-radiative processes. 

The above analysis holds when the overall system size 
is smaller than the radiation wavelength. This condi- 
tion for cooperative emission also allowed us to use the 
long-wave approximation for the Mie theory Green di- 
adic, where we disregarded the real part containing di- 
rect dipole-dipole interactions between the emitters. In 
absence of a NP, the latter can be considerably stronger 
than radiative coupling and may lead to a suppression 
of the photonic Dicke effect for large ensembles [15]. In 
contrast, in the presence of a NP and close to the SP res- 
onance, the self-energy matrix is dominated by its imag- 
inary part, while the dipole-dipole interactions lead to 
a sp read of dipole frequencies around average value loq 
[Ull^l. However, a such a disorder affects only sub- 
radiant modes by removing the degeneracy in their spec- 
tral positions, but has no significant effect on SR modes 
[l7l |. Therefore, the plasmonic Dicke effect can survive 
the dipole-dipole interactions even for large ensembles. 

Finally, the predicted plasmonic Dicke effect could be 
observed in experiments with controllable separation of 
emitters from the NP surface. In the recent experiment 
on Cy5 dyes linked with an Au NP ^ , a systematic study 
of distance dependence for the ensemble fluorescence was 
performed. Even though the number of emitters was not 
fixed, a fast saturation of quantum efficiency with in- 
creasing distance was observed, consistent with Eq. ([1]). 
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